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Background 

Until recently expression of housekeeping genes is presumed 
to remain constant in different cell types, across developmen- 
tal stages and various experimental conditions hence they are 
employed as internal controls to normalize gene expression 
and protein content data. Many studies conducted recently 
on plant and animal samples have however, reported variable 
expression of constitutive genes during cellular development 1 
and experimental conditions 25 . The reason for such variation 
may be that they are not exclusively involved in performing the 
basic function of the cell but they also contribute in other cel- 
lular processes 67 . The variation in their expression is therefore 
likely and it is advised to consider prior validation of the ex- 
pression stability of these reference genes, a routine step in all 
laboratories to accurately normalize the target gene(s) expres- 
sion. Normalization, a fundamental step is used to compare 
the target gene expression with the constitutively expressed 
housekeeping genes. Through this process the technical and 
biological variations that are introduced during sample prepa- 
ration or any other step during the experiment are controlled. 

In this report, the authors evaluated a panel of candidate 
reference genes for their expression stability in glioblastoma 
tissues using statistical tools namely Genorm 8 , Normfinder 9 
and equivalence test 10 . The hypothesis was to identify a set of 
housekeeping genes with minimal variability in glioblastoma 
tumors compared to non-neoplastic white matter. In addition, 
the significance of normalization by measuring the transcrip- 
tion levels of twelve target genes using different reference 
genes is demonstrated. 

Study design 

The study was intended to validate a panel of reference genes 
for their expression stability in glioblastoma (GBM) samples, 
the most common type of tumor in the central nervous sys- 
tem. These malignant gliomas are the most common and most 
aggressive type of glial tumors. Due to several complicating 
factors treatment of these patients is difficult. Thus, there is 
an urgent requirement to identify novel therapeutic targets 
for its therapy. In this context, many recent studies are con- 
ducted to explore high-throughput gene expression profiles us- 
ing real time PCR in order to identify novel genes differentially 
expressed in GBM. 



The authors first extracted mRNA from thirty glioblastoma tis- 
sues and nine non-neoplastic white matter tissue samples and 
subjected them to quantitative PCR. The 2-AACT equation was 
applied to calculate the relative expression of tumor samples 
and non-neoplastic brain tissues. To evaluate stably expressed 
reference genes in this experimental set up, the authors chose 
seven frequently used housekeeping genes (HKGs) from lit- 
erature. Candidate HKGs assessed were beta-actin (ACTB), 
glyceraldeyde-3-phosphate dehydrogenase (GAPDH), beta- 
glucuronidase (GUSB), hydroxymethylbilane synthase (HMBS), 
hypoxanthine guanine phosphoribosyl transferase 1 (HPRT1), 
TATA-box binding protein (TBP) and 18S ribosomal RNA (18S 
rRNA). The gene expression stability was assessed using pub- 
licly available softwares namely Genorm, Normfinder and 
equivalence test. Genorm program defines the gene stability 
as the average pairwise variation of a particular gene with all 
other control genes and ranks the genes according to their av- 
erage expression stability denoted by M. The gene with mini- 
mum M value is considered to be highly stable whereas the 
gene with highest M value is least stable and can be excluded. 
An alternative program, Normfinder, ranks the candidate ref- 
erence genes based on the combined estimates of both intra- 
and intergroup variations. The lower stability values indicate 
the more stable expressed candidate genes. The equivalence 
test estimates the significance of differences in the median 
expression values between tumor and normal tissues of each 
individual gene. 

For comparison of housekeeping transcription levels, the cycle 
threshold (CT) values were plotted directly, assuming similar 
threshold for all genes assessed. The seven housekeeping genes 
analyzed showed a wide expression range, with CT values be- 
tween 14 and 32. The CT is defined as the number of cycles 
needed for fluorescence to reach a specific threshold level of 
detection and is inversely correlated with the amount of RNA 
template present in the reaction. The CT values for all these 
genes showed normal distribution according to the Kolmogorov 
and Smirnov method, in both tumor and non-neoplastic sam- 
ples. These genes were clearly distributed into different expres- 
sion level categories. The extremely abundant 1 8S rRNA, which 
represents the bulk of total RNA in the cell, presented CT values 
below 18 cycles; genes coding for highly expressed mRNAs, 
such as ACTB and GAPDH, with majority of CT values between 
18 and 22 cycles; and the moderately expressed genes, GUSB, 
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HMBS, HPRT1 and TBP, showing CT values between 26 and 
30 cycles. Additionally, the authors found slightly higher levels 
in the abundance of GUSB and HPRT1 mRNAs than in HMBS 
and TBP mRNAs. 

The authors reported that amongst seven frequently used 
housekeeping genes TBP and HPRT1 were adequate references 
for glioblastoma gene expression analysis and that expression 
levels of ACTB, GAPDH, GUSB and 18S rRNA, differed signifi- 
cantly between tumor and normal tissues on the basis of the 
examination of raw CT values. This indicates that these genes 
are regulated and cannot be utilized for target gene normaliza- 
tions. It is recommended that this initial analysis is a crucial 
step in order to exclude highly unstable genes from further cal- 
culations. Based on data from a previously analyzed microar- 
ray, the authors selected 12 target genes signified as TG 1-12, 
candidates of being over-expressed in GBM, for quantitative 
RT-PCR experiments to investigate the correctness of overex- 
pression suggestion. The authors found five genes whose ex- 
pression levels significantly differ between tumor and normal 
samples. Surprisingly, two of them were decreased (TG4 and 
TG6 and three (TG1, TG8 and TG12) were increased in tumor 
tissue. Furthermore, Holliday Junction Recognizing Protein, a 
novel DNA repair protein over expressed in lung cancer, is found 
to extremely over-expressed in GBM samples. Evaluation of the 
expression levels of these genes utilizing different endogenous 
controls revealed that the normalization method applied might 
introduce errors in the estimation of relative quantities. Genes 
presenting expression levels which do not significantly differ 
between tumor and normal tissues can be considered either 
increased or decreased if unsuitable reference genes are ap- 
plied. Most importantly, genes showing significant differences 
in expression levels between tumor and normal tissues can be 
missed indicating the significance of prior validation of can- 
didate housekeeping genes for any experimental set up in a 
laboratory, especially when small differences are intended to 
be detected. 

Implications 

The authors found a fluctuated expression of Gapdh, beta- 
actin and 785 rRNA in glioblastoma samples which are fre- 
quently used as internal controls in majority of the laboratories. 
It is now demonstrated that expression of housekeeping 



genes, which are used as endogenous controls for normaliz- 
ing mRNA transcripts, may not be necessarily stable in all cells/ 
tissues under all conditions. A gene may show stable expres- 
sion in one condition whereas unstable in another. Therefore 
it is indispensible to pre-validate the expression stability of 
reference genes to accurately normalize the gene expression 
data. 
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